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•i-H 
Q^ Abstract 

Recently a model of intra- and interspecific competition between two species was proposed [Phys. Rev. 

E 87 (2013) 010101], in which the scarcer species (i.e., with smaller stationary population size) can be more 
\Q ■ 
l/~) ' resistant to extinction when it holds a competitive advantage. Here we verify this survival of the scarcer in 

G\ '. 

I/-, | space (SSS) phenomenon in models with spatial structure, both analytically and numerically. We find that 

the conditions for SSS, as obtained applying renormalization group analysis and Monte Carlo simulation to 

o ■ 

{•f. . a discrete-space model, differ significantly from those found in the spatially homogeneous case. 



Keywords: Demographic stochasticity, Doi-Peliti mapping, Renormalization group, Monte Carlo simula- 
tion. 

1 Introduction 

Recently, Gabel, Meerson, and Rcdncr [1. proposed a stochastic model of two-species competition in which, for 
stationary regimes such that the two species coexist, the scarcer species (i.e., with smaller population size), can, 
under certain conditions, be less susceptible to extinction than the more populous species. This phenomenon, 
dubbed survival of the scarcer (SS) in pQ, is surprising since conventional wisdom on population dynamics 
suggests that a smaller population is more susceptible to extinction due to demographic fluctuations. The 
authors of pQ use a variant of the Wentzel-Kramers-Brillouin- Jeffreys (WKBJ) approximation to obtain the 
tails of the quasi-stationary probability distribution of the process [21 E] ■ They show that a small asymmetry 
in interspecific competition can induce survival of the scarcer species. 



While the analysis of [TJ holds for well mixed (spatially uniform) populations, it is natural to inquire whether 
the same conclusions apply for populations distributed in space. In this paper we study a stochastic model 
similar to that of [JJ, but with organisms located on a d— dimensional lattice, and investigate the conditions 
required for survival of the scarcer in space (SSS). On reproduction, a daughter organism appears at the same 
site as the mother; competition only occurs between organisms at the same site. In addition to these reactions, 
organisms also diffuse on the lattice. Starting from the master equation, we build a representation involving 
creation and annihilation operators, taking advantage of the discrete nature of the stochastic process. This 
representation allows us to obtain a full functional integral formulation, which can be regarded as a mapping 
of the stochastic process to a field theory. This now standard procedure is known as the Doi-Peliti mapping 
[H[SJH]. We consider the population dynamics in a critical situation, i.e., populations close to extinction. The 
dynamic renormalization group (DRG) is used to study the nonequilibrium critical dynamics of the population 
dynamics, specifically, to determine how the model parameters transform under changes in length and time 
scales. The model is also studied via Monte Carlo simulation, which again confirms the existence of SSS when 
the less populous species is more competitive. 

The remainder of this paper is organized as follows. In section [2] we describe the model and write the 
effective action associated with the stochastic process. A brief analysis of mean-field theory, valid in dimensions 
greater than the critical dimension is performed. In section [3] we use some known results on multispecies 
directed percolation to obtain the renormalization group flow in parameter space. This analysis, combined 
with numerical simulations of the associated stochastic partial differential equations, allow us to determine the 
regions in parameter space in which SSS occurs. Results of Monte Carlo simulations are presented in Section 
|U Section \E\ closes the paper with our conclusions. 

2 Model 

The model proposed in [I] may be described, with some minor modifications in the rates as specified below, by 
the following set of reactions: 

A -3> A + A B A b + B 

A + A SL (A) B + B & (B) (1) 

A + B -^ B A + B -^ A 

The reactions in the first line correspond to reproduction of species A and B. Those in the second line 
describe intraspecific competition processes which may be of mutual annihilation (0) or individual death (A 
I B), while the third line represents interspecific competition, a (/?), a' ((3') and C (£) are the rates of the 
reactions associated with species A (B). 

The relation between the rates in equation ([TJ and those of the original model [1] is (in the case of mutual 
annihilation): a <+ 1, /3 <->• g, a' ■<-> 1/K, 0' <+ 1/K, ( <+ e/K, and £ O Se/K. 5 < 1 implies an asymmetry in 
competition between species that favors B at expense of A. In this case species B is more competitive than A; 
the opposite occurs if 5 > 1. 

The authors of Ref. [TJ constructed a phase diagram in the g-S plane (shown schematically in Figure [TJ. 
The phase diagram includes results of both mean-field theory (MFT) and analysis of the master equation for a 
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Figure 1: Phase diagram of well mixed system in the g-5 plane, as furnished by the analysis of Ref. pQ. Dashed 
blue curve: g m f = (1 + <5e)/(l + e). Full red curve: g sw = (1 + m<5e)/(I + me). 

stochastic population model in a well mixed system, i.e., without spatial structure. In MFT, the populations of 
the two species are equal when g m f(S) = (I + <5e)/(I + e). Analysis of the master equation yields equal extinction 
probabilities when g sw (S) = (1 + m<5e)/(I + me) with m = [(ln2) _1 — ll . These conditions are plotted in 
figure [TJ they intersect at the point (g, S) = (1, 1). The phase diagram includes two "normal" regions, in which 
the more populous species is less likely to go extinct, and two "anomalous" regions, in which the scarcer species 
is more likely to survive. Thus in region / (where <5 < I and B is more competitive in the sense defined above), 
species A is more numerous in the quasi-stationary state and simultaneously more susceptible to extinction, 
while in region // (where S > 1 and A is more competitive), species B is more numerous, yet more likely to 
go extinct first. Regions I and // are anomalous since they correspond to an inversion of the usual dictum 
that susceptibility to extinction grows with diminishing population size. Our goal in this work is to determine 
whether a stochastic model of two-species competition with spatial structure, in which organisms diffuse in 
d— dimensional space, exhibits a similar phenomenon. In the following subsection we develop a continuum 
description for this spatial stochastic process. 

2.1 Effective action 

We consider a stochastic implementation of the reactions of equation (TTJ) on a d-dimensional lattice, adding dif- 
fusion (nearest-neighbor hopping) of individuals of species A and B at rates Da and Db respectively. Following 
the standard Doi-Pcliti procedure mentioned in the Introduction, ignoring irrelevant terms (in the renormal- 
ization group sense) and with the so-called Doi shift .6^ already performed, we obtain the following effective 
action: 



S(<p, <f>, tl>,i/>) = / d d x / dt {<j>[d t + D A {a A - V 2 )]0 - a0 2 + jc/00 2 + a> 2 2 + CV>#} 

+ f d d x f dt {tP[dt + D B {a B - V 2 )]^ - (3^ 2 i> + j(3'^ 2 + /3'i> 2 ^ 2 + f #$} , (2) 



where a a = —ci/Da, cb = —ft/Ds and J = 1 (2) for individual death (mutual annihilation). Action S has 
four fields. The expected values of the fields <f> and tp (evaluated by performing functional integrals over the 
four fields, using the weight exp[— S]), represent the mean population densities of species A and B, respectively. 
The fields with overbars have no immediate physical interpretation, but are related to intrinsic fluctuations of 



the stochastic dynamics. 

Using the definitions of a a and <7b and the change of variable (f> — ¥ 9a4>, 4> ~ ► 
where 9a = 



M£- , we can write the effective action in the form: 



<f>,ip—} Ob^i and ip 



iV. 



with 



Dimensional analysis yields, 



S(^,0,V>,VO = J d d x J dt^[d t +D A (a A -V 2 M 
+ ^[d t +D B (<j B -V 2 )}^ 

+ h B H J $ + hA(f>ip4>\, 
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[/ia] = [hi 



„2-# 



l5Uj = [9b\ 



(3) 



(4) 



where [X] denotes the dimensions of X, and p has dimensions of momentum [BJ. The upper critical dimension 
is therefore d c = 4, as for single-species processes such as the contact process and directed percolation. 



2.2 Mean-field approximation 

For d > d c , mean-field analysis, which ignores fluctuations in <p and ?/>, yields the correct critical behavior. 

6d> 



Imposing the conditions |§ = 0, t§ = 0, |4 = and ~ = 0, we obtain the mean-field equations 



ff = D A V 2 + a^-g A <j> 2 - 

72,/, i O.I, _ „„„/,2 



(5) 



ff : D B V^ + ^-gB^-h A H- 

Let L be a typical length scale and define X = gA<t>/a, Y = gsip/ft, s — DAt/L 2 , and define a rescaled 
coordinate x' = x/L. Then equation ((5]) can be written in dimensionless form as 



(6) 



^ = DV 2 Y + -/(cY-cY 2 -bXY) = DV 2 Y + g(X,Y) 

with 7 = aL 2 /D A , a = (3h B /(ag B ): b = /3h A /(ag B ), c = g A (3 2 /(g B a 2 ), D = D B /3gA/(DAag B ), f(X,Y) = 
7 [X — X 2 — aXY) and g(X, Y) = 7 (cF — cF 2 — WfYJ . Coexistence is possible in the stationary state if 
the conditions c > b and a < 1 hold, implying /i^ < /3gA/a and /is < ags/P- These conditions depend 
monotonically on parameters a, /?, 3,4 and gs 1 analogous to the coexistence conditions found in pQ. The same 
is not true in spatial stochastic theory, as we shall see in the next section. 



3 Renormalization group (RG) flow 

Some years ago, Janssen analyzed a class of reactions of the kind defined in equation ([!]) [7 . This work 
considered multi-species reactions of the form 

Xi^2X t X t ^% Xi+Xj -+kXi+lXj, (7) 

where i and j are species indices and k, I are either zero or unity; this process is called multicolored directed 
percolation (MDP), with different colors referring to different species. The RG analysis of [7J shows that 
knowledge of the directed percolation fixed points is sufficient to determine the fixed points of the full MDP. 
As shown in that work, the parameter combinations gAJ Db = ua and gs/ Da = ub, related to intraspecific 
competition, flow under renormalization group transformations to the stable DP fixed point u* A = u* B = u* = 
2e/3 with e = 4 — d > Oo Therefore, regarding the renormalization of intraspecific competition parameters, 
inclusion of other species behaving as in DP does not alter the fixed point. 

Janssen's analysis shows that there are four fixed points for interspecific competition parameters Ha/ Da = 
va and hs / Db = vb, which from reactions (JTJ) are related to £ and £. The first two are (y A , v B ) = (0, 0), which 
is unstable, and 

, * M (Da + Db2 D a + Db 2\ 

which is a hyperbolic fixed point if Da = ObO For Da ^ Db, it was conjectured [7] that the stability of the 
fixed points remains the same despite small changes in the renormalization group flow diagram topology. The 
other two fixed points for Da = Db are stable; they are given by (v A ,v B ) = (0, 2u*) and {v A ,v B ) — (2u*,0). 
To summarize, for Da = Db = Dq, the interspecies competition parameters flow to one of the following four 
d— dependent fixed points (denoted as 0,F,G, and H below), 



(va,Vb) 



tft 9bH 



D Q ' D 
{O:(0,0), F :(«*,«*), 

G:(2«*,0), H:(0,2w*)} (9) 



Figure [5] shows the renormalization group flow diagram. The unstable fixed point O corresponds to complete 
decoupling between the two species, and the hyperbolic point F to symmetric coupling. In the symmetric 
subspace, any nonzero initial value of v = va — Vb, flows to F. For vb > va (S > I), (va,vb) flows to G, so 
that species A is effectively unable to compete with B. Similarly, if va > vb (8 < 1), the flow attains point H. 
In the following subsection we discuss the stationary states associated with these fixed points. 



1 Not to be confused with parameter e used in equation [T] 
2 In the nomenclature of [7|, if species are of the same flavor. 




Figure 2: RG flow of the interspecific competition parameters vb and va- Only the first quadrant is relevant 
due to the requirement of positive rates. There is an unstable fixed point O at (va,vb) — (0,0) (black dot); 
a hyperbolic fixed point F at (va,Vb) = (u*,u*) (red dot), and two stable fixed points G and H, located at 
(va,vb) — (0,2m*) and (va,vb) = (2w*,0) respectively (blue dots). The separatrix va = vb is the boundary 
between the basins of attraction of G and H. 

3.1 Steady state close to critical points 

Population dynamics in the critical regime is governed by a pair of coupled stochastic partial differential equa- 
tions (SPDE), which are readily deduced from the action of equation ([3]) 8 : 



dt 



.DqV cj) + a<f>- g A <t> - h B <H + Vi 



(10) 



and 



-£ = D VV + M - 9B^ 2 - h A ilxt> + V2 
at 

where the noise terms ?yi(x,i) and 772 (x,i) satisfy (771 (x, £)) = (rj2(~x.,t)) — 0, and 



(11) 



(r7i(x,i) m (x,t))=2 ffA ^(x,t)^(x-x')<5(t-t') ) 



(12) 



(r? 2 (x, t)r? 2 (x, i)) = 2 5B ^(x, t)5 d {* - x')8(t - t'). 



(13) 



Note that these multiplicative noise terms depend on the square root of their respective fields, and were obtained 
directly from the action, without any additional hypotheses. 

3.1.1 SPDE in the vicinity of point F 

In the symmetric subspace, (5=1, the RG flow is to the fixed point F, and parameters qa/Dq, 5b/-Do, Ii-a/Do, 
and hs/Do take the associated values. Therefore the SPDEs in (fT0| and (fTTj) can be written as [9]: 



-£ = A)V 2 + (jy - D u*(b 2 - D u*<H; + r h , 
at 



(14) 



^=D V 2 ^ + g^-D u*^ 2 -D u*^ + r l2 (15) 

at 

where we set a = 1 and /3 = g. With the rescaling (f> — >• <j>/(Dqu*), ip — ► ^/(Dqu*) and x — ► (l/Do) 1 ^ 2 ^, we can 
write equations (fT4| and ([T5j) in the form: 

^ = V 2 + 0-^-^ + 7 ?1 , (16) 

|=V^ + #-^ 2 -# + r, 2 (17) 

with rescaled noise 

( ( , 1 (x,t)i, 1 (x,())=2(fl ) 2 -M Il ')^(x,^(x-x>(t-O (18) 

(»fc(x,t)»fe(x,t)) =2( J D ) 2 -M U *) 2 ^(x,i)^(x-x')<5(t-O- (19) 

3 

For d < 4, the rescaled noise intensity <r grows with diffusion rate Dq; for d = 1 we have a 2 = 2D^u* 2 with 



u* 2 



2. Figures 3(a) and |3(b)j show results of numerical simulations of equations (J16I) and (IT7l) in d = 1 with 



g = 1. They show, for a one-dimensional lattice of L = 128 sites in the interval (—1, 1), the population densities 
averaged over 1000 Monte Carlo realizations for two different values of the noise intensity. The time interval 
T is T = 400, partitioned into N = 40000 steps so that A* = T/N = 400/40000 = 0.01. Initial conditions 
are (/>(x,0) = 0.4 and ip{x,0) = 0.6. For a — 0, population densities are almost time-independent as shown in 



Figure 3(a) This is not the case for larger noise values. In this case, the two populations tend to the same 
value (see figure [3(b)] ) . Although Figs. 3(a) and |3(b)| represent averages over many realizations, we have verified 
coexistence in individual runs. 

These simulations were performed using standard integration techniques for the Langevin equation with 
the XMDS2 software 10]. Since we are interested in showing only the initial temporal trends of the two 
population densities in the vicinity of fixed points, the use of standard techniques of integrating the Langevin 
equations is sufficient to reveal the qualitative nature of the solutions. Near a phase transition to an absorbing 
state, one of the population densities tends to zero and the standard numerical integration scheme fails. This 
standard algorithm is not suitable for extracting the more accurate results associated with critical exponents or 
asymptotic decays. In this case we would have to use more sophisticated algorithms such as those proposed in 

punned]. 

3.1.2 SPDE in the vicinity of a stable fixed point 

Now consider the case 5^0. With S > 1, the parameters flow to point G and Eqs. (TIT)]) and (TTTj) become 



^ = A)V 2 + 0-^ + 77!, (20) 

^^D V 2 ^ + g 1 p-u* 1 p 2 -2u*^ + 7 l2l (21) 

at 

where we put a = 1 and /3 = g. With a rescaling similar to that used above, we have 
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(a) Symmetric case with a = 0.0 and g = 1. Initial popu- (b) Symmetric case with a = 0.03 and g = 1. Initial popu- 
lations are </>(x, 0) = 0.4 and f/>(a:, 0) = 0.6. lations are tf>(x, 0) = 0.4 and ^(x, 0) = 0.6. 

Figure 3: Numerical simulations symmetric SPDE. 
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at 



(22) 

(23) 



with 771 and 772 satisfying (fT8|) and (|19p . respectively. For (5 < 1 the situation is completely analogous to the 
case 8 > 1, and the resulting equations are as above with the exchange of and ?/> and changing the position of 
the factor g accordingly. 

If we neglect diffusion and noise terms, we have a system of ordinary differential equations having a fixed 
point associated with the stable coexistence state given by ((j>*,ip*) = (1,5 — 2)|f| Therefore, the condition for 
coexistence is g > 2. If g < 2, only species A survives for 8 > 1. Figure (0} shows the result of a simulation 
with g = 2.2 in which the equations are integrated including both diffusion and noise. Numerical experiments 
indicate that the effect of these terms is to increase the value of g ^ 2 for coexistence. Higher noise values imply 
higher thresholds g for coexistence. 
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Figure 4: Asymmetric case (8 = 1.5) with a — 0.03 and g — 2.2. Initial populations for A and B species are 
4>(x,0) = 0.1 and tj){x,Q) = 0.9, respectively. 

For 8 < 1 similar reasoning applies. The results are summarized in the phase diagram of Fig. [SJ There are 
four stationary phases: a pure-A phase for 8 > 1 and g < 2, a symmetric pure-i? phase for 8 < 1 and g < 1/2, 
and two (disjoint) coexistence phases. The latter arise when the less competitive species proliferates sufficiently 
faster than the more competitive one. 



'For the case S < 1, fixed points are (<f>*,ip*) = (1 — 2g,g), and coexistence condition is < g < 1/2. 
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Figure 5: Stationary phase diagram in the g-S plane for stochastic model with spatial structure. The reference 
points have coordinates: b = (1,0), p = (0, 1/2), r = (3,2). For 6 > 1, when g < 2 only species A is present, 
while for g > 2 there is coexistence. Similarly, for 5 < 1, when 5 > 1/2 only species B is present, and for g < 1/2 
there is coexistence. The diagonal line connecting points p and r is the equal-population criterion furnished by 
mean-field theory, g m f(S) = (1 + e<5)/(l + e), with e = 1. In regions / and II, the more populous species (A 
or B, respectively), as predicted by mean-field theory, in fact goes extinct. The thick, dashed, blue and black 
curves represent g m f and g sw for the value of e that maximizes the area between them, i.e., e = l/^/m w 0.66. 

3.1.3 Nature of phase transitions 

From the equations that omit the diffusion and noise terms, one may infer the nature of the phase transitions in 
the diagram of Fig. ([3]). From Table [TJ which shows the fixed points of the equations for different values of 5, we 
see that in general the values in the final column do not match for 6^1. This fact indicates that the vertical blue 

Table 1: Simplified equations and their fixed points for different values of 5 



5 


Equation 


Fixed Point 


= 1 


<f> = <f> — (j) 2 — (fitp and ip = gip — ip 2 — <pip 


(P,p) = (0,g) 


>1 


cb = <p — 4> 2 an d ip = gip — ip 2 — 2<j)i(j 


(4>*,P) = (l,2-g) 


< 1 


cb — g4> — cj) 2 and ip — ip — t/j 2 — 2(bip 


{<P,<P) = {g,l-2g) 



line in the diagram ([5]), i.e., the line 5 = 1, is a line of discontinuous phase transitions. Similarly, examining the 
expressions in the final column of Table [TJ we infer the nature of phase transitions along the horizontal lines at 
g = 1/2 and g = 2. In regions of coexistence, the densities of the minority species grow continuously from zero; 
thus these transitions are continuous. One should of course recognize that the predictions for the phase diagram 
are essentially qualitative; quantitative predictions for nonuniversal properties such as phase boundaries are 
not possible once irrelevant terms have been discarded. Regarding the nature of the transitions, we expect the 
continuous ones, as furnished by this mean-field-like analysis, to in fact belong to the DP universality class 
[7] . The line of discontinuous transitions contains a portion (between levels p and r) that is between absorbing 
subspaces (only on species present in each phase), and other portions that separate (from the viewpoint of the 
species undergoing extinction) an active and an absorbing phase. Discontinuous transitions of this nature are 
not expected to occur in one dimension [14] , 

According to mean-field theory, a point in region iT corresponds to ps > Pa- If we increase the competi- 
tiveness of species A, increasing S while maintaining g fixed, there will be a point beyond which species B no 



longer has the greater population density. Improvement of the competitiveness of species A will make it the 
more populous species if organisms are well mixed. In the case of a spatially structured population, by contrast, 
any competitive advantage of species A, no matter how small, is sufficient to dominate the B population (if 
the latter does not proliferate very rapidly) when the population densities are very low, as is the case near 
criticality. Thus in a situation in which mean-field theory predicts species B to be the majority, we can have 
species A exclusively. This can be seen as a strong version of the survival of the scarcer phenomenon. 

To compare the predictions of the mean-field and spatially structured descriptions, a pair of dashed lines 
are plotted in Fig. [5l representing equations for g m f and g sw as defined previously, for parameters such that 
the area between them is maximum, thereby maximizing the region in which SS occurs. [The maximum area 
is obtained using e = 1/y/m with m = [(ln2) — 1] , in Eq. (TTJ]. Given the much larger area of the region 
exhibiting SS in the spatial model, compared with that in the mean-field analysis (i.e., the area between the two 
dashed lines for 1 < 6 < 3,) we see that the SS phenomenon can be greatly intensified when spatial structure 
is included. Analogous reasoning applies for region / (green). (For the parameter values used, however, this 
reasoning does not apply for S > 3, since in this case the area between the dashed lines can be arbitrarily large. 
Also we no longer have a region analogous to region II with species A only.) In this way we see that the phase 
diagram for the spatial stochastic model is rather different from that found in [I]. Moreover, under certain 
conditions, the survival of the scarcer phenomenon is strengthened. 

4 Simulations 

Since the preceding analysis involves approximations whose reliability is difficult to assess, we perform simula- 
tions of a simple lattice model to verify SSS. We consider the following spatial stochastic process, defined on a 
lattice of L d sites. Each site i is characterized by nonnegative occupation numbers Oj and bi for the two species. 
The organisms (hereafter "particles") of the two species evolve according to the following rates: 

• Particles of either species hop to a neighboring site at rate D. 

• Particles of species A(B) reproduce at rate Xa(^b)- 

• At sites with two or more A particles, mutual annihilation occurs at rate aAO,i(a.i — 1), and similarly for 
B particles, at rate asbi(bi — 1). 

• At sites having both A and B particles, the competitive reaction B — ► occurs at rate ^A^ibi and the 
reaction A — ► occurs at rate C,BO-%bi- 

The model is implemented in the following manner. In each time step, of duration At, each process is 
realized, at all sites, in the sequence: hopping, reproduction, annihilation, competition. 

In the hopping substep, each site is visited, and the probability of an A particle hopping to the right (under 
periodic boundaries) is set to ph = Da, At/ (2d). A random number z, uniform on [0,1) is generated, and if 
z < ph the site is marked to transfer a particle to the right. Once all sites have been visited, the transfers are 
realized. The same procedure is applied, in parallel, for B particles hopping to the right. Subsequently, hopping 
in the other 2d — 1 directions is realized in the same manner. 



10 



In the reproduction substep, the A particle reproduction probability at each site i is taken as p c = A^a^Ai. 
A random number z is generated and a new A particle created at this site if z < p c . An analogous procedure 
is applied for reproduction of B particles. 

In the annihilation substep, a probability p a = aA<ii(ai~ 1) At is defined at each site, and mutual annihilation 
(a,i — > a,i — 2) occurs if a random number is < p a . Again, an analogous procedure is applied for annihilation of 
B particles. 

Finally, at sites harboring both A and B particles, probabilities pa = QA^ibiAt and ps — CBO-ibiAt are 
defined. The process B — > occurs if a random number z < pa, while the complementary process A — > occurs 
if PA < z < pa + Pb- Note that at most one of the processes (A — ► and B — \ 0) can occur at a given site in 
this substep. 

The time step At is chosen to render the reaction probabilities relatively small. To do this, before each 
step we scan the lattice and determine the maximum, over all sites, of Oj, bi, and a^. With this information 
we can determine the maximum reaction rate over all sites and all reactions. Then At is taken such that the 
maximum reaction probability (again, over sites and reactions) be 1/5. In this way, the probability of multiple 
reactions (e.g., two A particles hopping from i to i + 1, etc.) is at most 1/25, and can be neglected to good 
approximation, particularly in the regime in which occupation numbers are typically small. 

We simulate of the process on a ring (d = 1) using quasistationary (QS) simulations, intended to sample 
the QS probability distribution [1 51 116] . In these studies the system is initialized with one A and one B 
particle at each site. When either species goes extinct (an absorbing subspace for the process), the simulation 
is reinitialized with one of the active configurations (having nonzero populations for both species) saved during 
the run. Following a brief transient, the particle densities fluctuate around steady values. We search for 
parameter values such that species A is more numerous, despite being much less competitive than species B 
(i.e., Cb » CO- 

Given the large parameter space, certain rates are kept fixed in the study: We set D — ola = &b = 1/4. 
We use QS simulations as well as spreading simulations [T7] to estimate the critical value of A for single-species 
survival. In spreading simulations the initial configuration is that of a single site with one particle, and all 
other sites empty. We search for the value of A associated with a power-law behavior of the survival probability 
P(t) and the mean population size n{t). These studies yield A c = 1.267(1) as the critical point for survival of 
a single species, i.e., without competition. (Here and in the following, figures in parentheses denote statistical 
uncertainties in the final digit.) The phase transition is clearly continuous; details on scaling behavior will be 
reported elsewhere. 

In the two-species studies we set A^ = 1.6 and (a — 0.005, so that species A is well above criticality but 
weakly competitive. We then vary As and £b, monitoring the QS population densities pA and ps of the two 
species, as well as their QS lifetimes, ta and tb- The latter are estimated by counting the number of times, in 
a long QS simulation, that one or the other species goes extinct. Survival of the scarcer is then characterized 
by the conditions pa > Pb and ta < tb- For the parameter values studied, we observe SSS for Xb ~ A c and 
Cb ^ Ct- Figure H2 illustrates the variation of the population densities and of the ratio ta/tb as Cb is varied 
for fixed Xb = 1.25, on a ring of 200 sites. In this case, pa > Pb throughout the range of interest; the scarcer 
species, B, survives longer than A for Cb > 0.25. Figure [7] shows a sample evolution, for Cb — 2 and other 
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parameters as in Fig. [6l illustrating that sites bearing both species are quite rare, occurring at the frontiers 
between regions dominated by a single species. 
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Figure 6: Quasistationary population densities pa (blue) and ps (red), and lifetime ratio ta/tb (black) versus 
Cs ; for parameters as specified in text. 

Figure \E\ shows the result of a systematic search in the Cb~Xb plane on a ring of L = 200 sites. These data 
represent averages of 200 realizations, each lasting 10 5 time units. SSS is observed in the region bounded by 
the two curves; the lower curve corresponds to ta = tb while on the upper we have pa — Pb- Thus SSS is 
found in a small but significant region of parameter space. Although we use a rather small system to facilitate 
the search, SSS is not restricted to this system size. For L = 800 for example, the range of Xb values admitting 
SSS is more restricted (from about 1.24 to 1.27) but at the same time the effect can be more dramatic. For 
Xb = 1-27 and (b = 0.2, for example, we find pa = 0.74, ps = 0.49, and tb ~ lOOr^. Thus SSS is observed in 
simulations of the lattice model. Studies using A^ = 1.35 and £a — 0.01 yield a qualitatively similar diagram, 
leading us to conjecture that the form of the region exhibiting SSS in one dimension is generically that shown 
in Fig. E 

In the inset of Fig. [5] the data are plotted in the S-Xb plane, for comparison with the theoretical predictions 
shown in Fig. [5] The region exhibiting SSS is rather different from the predictions of both MFT and stochastic 
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Figure 7: A sample evolution on a ring of 200 sites, of duration 4000 time units (with time increasing downward). 
Blue: sites with c^ > and bi = 0; red: sites with a, = and bi > 0; black: sites with both species present. 

field theory. In particular the lower boundary in the S-Xb plane is not horizontal and it is unclear whether 
the SSS region extends to 5 = 1. (As Xb is increased, SSS is observed in an ever more limited range of Cb 
values, making numerical work quite time-consuming.) Here it is important to recall that irrelevant terms (in 
the renormalization group sense) are discarded in the theoretical analysis, so that predictions for the phase 
diagram are only qualitative. More detailed simulations, including other parameter sets, larger systems, and 
simulations on the square lattice, are planned for future work. 

5 Conclusion 

We study a spatial stochastic model of two-species competition, inspired by the spatially uniform system ana- 
lyzed in Ref. [I] . Based on Janssen's results of for multicolored directed percolation [7] , we find that the flow of 
parameters under the renormalization group is such that a small competitive advantage may be amplified to the 
point of excluding the less competitive species, unless the latter reproduces rapidly. Thus the less competitive 
population can go extinct, regardless of its size as given by mean-field theory. These effects tend to be stronger, 
the smaller the dimensionality d of the system. Our result shows that the survival of the scarcer phenomenon 
observed in [I] persists, in some cases in a stronger form, in the presence of spatial structure. These differences 
are evident when comparing Figs. ((TJ) and ([5]). 

Monte Carlo simulations of a one-dimensional lattice model verify the existence of SSS when the more 
competitive species has a reproduction rate near the critical value, while that of the less competitive one is 
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Figure 8: Simulations in one dimension: SSS is observed in the region bounded by the two curves; see text for 
parameters. Inset: the same data plotted using S = Ca/Cb instead of Cb on the horizontal axis. 

above criticality. 
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